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' We present an analytical study of one-dimensional semiconductor superlattices in external electric 

, fields, which may be time-dependent. A number of general results for the (quasi)energies and 

eigenstates are derived. An equation of motion for the density matrix is obtained for a two-band 
model, and the properties of the solutions are analyzed. An expression for the current is obtained. 
'. Finally, Zener-tunneling in a two-band tight-binding model is considered. The present work gives the 

background and an extension of the theoretical framework underlying our recent Letter [J. Rotvig 
et at, Phys. Rev. Lett. 74, 1831 (1995)], where a set of numerical simulations were presented. 

73.20.Dx,73.40.Gk,73.50.Fq 

> 

O ' I. INTRODUCTION 

o ■ 

. Studies of Bloch electrons under the influence of external electric fields attract at present intensive theoretical 
\^ ' attention (some recent papers can be found in Refs. ]l[-p^); this is a natural consequence of recent experimental 
^\ ' advances, which include the observation of Bloch oscillations , and studies of photon-assisted transport [p"5|-p7[ . 

• Thus the classic predictions of Bloch |Q and Zener have finally been verified, and the arena is open for new 
Ci I investigations and ideas. The physics of superlattices in external fields is extremely rich due to the large number 
^ . of parameters that can be controlled quite freely. Many of the physical properties are sensitive functions of these 
parameters, and hence slight adjustments in their values allow one to move between different physical regimes, both 
experimentally and theoretically. Examples of such parameters are the miniband structure of the superlattice (which 
can be controlled by varying the composition and/or thickness of the layers comprising the superlattice to meet the 
^ , requirements of the particular investigation) or the intensities and frequencies of the external fields. This flexibility 
allows one to study many different physical phenomena. The papers quoted above have addressed such varying physical 
phenomena as (i) dynamical localization or band-collapse (originally discussed in |2^) either in dc or ac electric fields, 
(ii) interplay of field-induced localization and Anderson localization due to disorder, (iii) chaotic motion of charge 
?H ' carriers, or (iv) full-scale numerical integration of semiconductor Bloch equations, which allow one to study interaction 
effects, such as exciton dynamics. 

Much of the interest has been caused by the need of understanding the interplay of Bloch oscillations {i.e., coherent 
time-periodic motion of charge carriers in one band), Zener tunneling between Bloch bands, and interaction effects. We 
studied recently |pl]] a two-band tight-binding model, and found via a direct numerical solution of the density-matrix 
equation-of- motion a rather complicated temporal behavior for the diagonal (in miniband index) component. This 
structure could, however, be interpreted in terms of intuitive physical arguments. Our numerical study was facilitated 
by a substantial amount of analytic background work, and the purpose of this paper is to provide a full account 
of the formal developments. In subsequent sections we shall derive a number of formal results, which we believe 
to be new, or generalizations of previously known results, and we give mathematical proofs of several statements 
made in our Letter. We also give a detailed discussion of the density-matrix equation-of-motion, upon which our 
previously reported numerical work was based, and which can be used as a starting point for future studies, for 
example consideration of relaxation or loss of phase coherence due to scattering. 

The paper is organized as follows. In Section II we derive the energy spectrum for a A^-band superlattice, in 
static and time-dependent fields, and apply the general results to the tight-binding two-band superlattice of Ref. pl| . 
Section III is devoted to the equation-of-motion analysis of the density matrix, and a derivation of an expression for 
the current, from which semiclassical intraband and quantum-mechanical tunneling contributions can be identified. 
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We also give some analytic results for Zener tunneling, which we used in the interpretation of our numerical results 
of Ref. The final section gives our conclusions. 



II. ENERGY SPECTRA 



In this section we consider the energy spectrum of an electron in a superlattice in the presence of a static or time- 
periodic electric field, which causes coupling between the bands. Our formulation allows us to obtain a simple proof of 
the existence of the Wannier Stark (WS) ladders in an arbitrary superlattice with band coupling, when a static field is 
applied. The number of interpenetrating WS ladders is equal to the number of zero-field energy bands in the model. 
This was first proven by Emin and Hart [|l2| by splitting the electric potential into two parts with a sawtooth and a 
staircase shape; our proof does not rely on this particular construction. In the final part of this section we introduce 
the tight-binding model of a superlattice containing two minibands used in the simulations of | |2l[ |, and establish a 
connection to the general results derived in the beginning of this section. 



A. Electronic motion in a superlattice 

We consider a semiconductor superlattice with a growth direction parallel to the x-axis. The lattice period is 
denoted by d and is repeated Np = L/d times, where L is the length of the superlattice. In the effective mass 
approximation |22| the wave function satisfies the usual Schrodinger equation with the scalar-potential Hamiltonian 
given by 

H4,{x, i) = 1^ + Vc{x) + eE{t)x , (1) 
Am 

where m is the effective mass, Vc{x) the periodic potential with period d, and E(t) (£'(t),0, 0) is the electric field. 
For E{t) = 0, the eigenstate problem Htp{x)4'{x) — eip{x) is solved by the usual Bloch states (BS) 

(f>nK{x) = -^e'^^'UnKix) (2) 

VL 

with eigenvalues en{K) associated with the miniband indexed by n. In a finite electrical field we follow and 
introduce accelerated Bloch states (ABS) when calculating the energy spectrum. We assume that the field is applied 
a.t t — and define the vector potential A{t) = {A{t), 0,0) by A{t) = — E{t') dt' . Then the semiclassical time 
evolvement of the crystal momentum k(t) = (fc(i), 0, 0) is given by 

k{t) = K ~ (e/h) [ E{t')dt' = K + eA{t)/h;. (3) 



At this point it is useful to introduce the Hamiltonian Ha in which the electrical field is represented by the vector 
potential instead of the scalar potential 

HAx,t)=^J^^±^ + V.ix). (4) 

The ABS, defined as 

1pnK{x,t) = -^e'^''Unk{t){x), (5) 

are instantaneous eigenstates to the time-dependent Hamiltonian Ha{x, t), 

Ha{x, t)%l)nK{x, t) = en[k{t)]lj)nK{x, t) . (6) 

The set of ABS forms an orthonormal basis. The so-called Houston function 4'nk{t) {x) can be expressed in the 
ABS basis as 

^ ^e^Ht)^u^f,^,j{x) . (7) 
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The Houston function is thus an ordinary BS corresponding to the wavevector k(t). 
We shall write solutions to the Schrodinger equation, 



H4x,mx,t)=zn^^l^ , (8) 



in the form ^{x, t) — e ' u{x, t), where u{x, t) is spatially periodic with period L. We thus have 

H^{x, t)u{x, t) = eu{x, t) + ■ (9) 

For a static electric field the function u can be chosen to be independent of time and e is the energy. If the field 
is periodic in time with a period Tac^ Floquet's theorem states that u will have the same periodicity. The Floquet 
solutions define the quasienergy e which is a generalization of the energy. 
Now let us expand ipix.t) in the ABS according to 

i^(x, t)^Y. CnK{t)e''^^'^^/''i;nK{x, t) . (10) 

nK 

We then insert (|l^) into (||) and use 

iTi — = in—j=-V kUnk(t)[x)k , (11) 

where we defined k = dk/dt. Since Vfcit„fc is a periodic function with period d, it may be expanded in terms of the 
functions Un'k according to 

n' 

where the expansion coefficients i?„'„ are the matrix elements 



Rn'n{k) = -I da;u*,fc(x)VfeU„fc(a;) . (13) 

" ■l-d/2 



Combining Eqs.( |ll[ - |l2| ) with k = —eE/h results in the following coupled equations for the expansion coefficients 

■f^ 9CnK{t) ^ ^^^m^^CnKit) " Fit) Rnn' [kit)]C„' K (t) , (14) 

n' 

where F{t) — —eEit). Equations similar to Eq.(|l^) were derived by Krieger and lafrate |^^. Our formulation differs 
from theirs, in that they use a gauge transformation in order to eliminate the spatially non-periodic term eE{t)x. 
We avoid the gauge transformation and thereby keep the notion of energy. Instead we introduce a common factor 
exp[ieA{t)x/fi\ in the ABS-expansion of the wavefunction. As seen from Eq.(pO|) this is the same as using an ordinary 
BS-expansion of the wavefunction: 

i^{x,t) = ^C„if(t)0„fc(t)(a;) 

nK 

*/^'^c„K(t)0„K(x) , (15) 

nK 

where CnK{t) = Cn^K-eA{t)/n{t)e'''^*^^ ■ The function u is found from Eq.(|5|) 

u{x,t) = ^CnK{t)(f>nKix) ■ (16) 
nK 

The set of coupled equations ( p^ for the coefficients CnK (t) will now be shown to lead to the existence of WS-ladders 
in static electric fields. 
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1. Static fields 



In the following we consider a model periodic potential corresponding to a finite number N of energy bands and 
demonstrate the existence of WS-ladders in static fields. For a static field, F(t) — F{> 0), the c„x's in Eq.(|l^) are 
independent of time. From this condition and the periodicity of the c„if 's in K it follows that 

CnK{t + Tb) = CnK{t)e-'''^-'f' , (17) 

where Tb = 2n/u!B and ljb = Fd/Ti. For each K in the Brillouin zone (BZ) there are linearly independent solutions 
Cifj(t) \ j = 1, ...,N to the system Eg. (p^) . From Eq.(p^) we see that the energy for a specific solution CKj{t) is 
only determined modulo Tllub- However, the substitution e ^ e + phuJB , P — integer, yields CKj{t) c/^^ (t)e'^"^^*. 
Only one of the corresponding solutions c^^it), p = integer, can therefore be time-independent. Let us denote the 
corresponding energy by exj ■ The system of equations (|lj) has periodic coefficients with period Tb ■ By applying the 
A''-dimensional Floquet theorem we know that the solutions are of the form CKj{t) = e"^^J*P5fj(t) , where lukj is 
only determined modulo lub and P0f^.(i + Tb) = P?f,(t). From Eq.(0) the energies of the Floquet solution C/^j(t) 
are = —htOKj + phujB , P — integer. Let p{K,j) be given by e^j^'"''' = exj- Then the energy spectrum is 

EKj = -hLUKj+p{K,j)hLJB ; j = l,...,iV ; KeBZ . (18) 

We can obtain more information by eliminating the K dependence in Eq . ([l^) : the semiclassical motion of k{t) is 
monotonic and can therefore be inverted to yield t{k) — h{k — K)/F. By use of the definition Dn{k) — c„fce^"'^/^ we 
obtain CnK{t) = £'„[fc(i)]e"^^/^. Eq.(|l^) is then reduced to a system of equations which is independent of the initial 
momentum 

iF^^^^^eJk)D^{k)~FY^R^,Ak)Dn'{k) . (19) 

n' 

A direct calculation shows that 

Dn{k + = D„(fc)e-"^«/'' . (20) 

As before the energy is only determined modulo Tilob- If for a definite solution T)j{k) to Eq.([l9|) we make the 
replacement e ^ e + phuiB , P — integer, then ckj — > ckjc'^^p'^. Therefore the displaced energy corresponds to a 
spatial translation of the solution by —pd, i.e., the eigenstates in terms of the Dj(fc)'s are only determined within a 
direct lattice displacement. The system Eq.([T9|) has periodic coefficients with period 2Tr/d. Appealing again to the 
A^-dimensional Floquet theorem we conclude that the solutions are of the form Dj(fc) = e*'^''j Pj(fc) ; j — 1, N . The 
Floquet functions P](fc) are periodic with period 2TT/d. By calculating Dj(fc + 2TT/d) and using Eq.(po|) we therefore 
obtain the energy spectrum 

'^jp = -F'^j +P^^B ; i = 1, N ; p integer . (21) 

We see that the spectrum can be visualized as N WS energy ladders. The relative positions of these are expressed in 
terms of the Floquet coefficients rj. This completes our general proof of the existence of Wannier Stark (WS) ladders 
in a static electric field. 

It should be noted that the lokj^s and r^-'s are connected. This can be demonstrated as follows. Floquet solutions 
for D(fc) can also be found by using the Floquet solutions for Cxit). For an arbitrary choice of K we put Dj(fc) = 
e-''<^K/FCKj{k) = e'^^K.k/ F ^~t(e+hu,K,)K/ F-pQ^_ij,y rpj^g energy spectrum is then e^^- -Tiujk] +pfiojB, P integer. 
Uniqueness of the spectrum implies that the energies {fuvKj \ j — i,---,N} are equal to the WS ladder positions 
{Fvj \ j — 1, ...,N} except for individual multiples of Hlob- 



2. Time-periodic fields 

In this section we consider time-periodic solutions u{x,t) to Eq.(^). We can derive general results concerning the 
quasienergy spectrum using the properties of the C^f-coefficients. As we shall see, the interplay between the spatial 
periodicity and the temporal periodicity can have significant consequences for the quasienergy spectrum. 

Specifically, let us consider a time-dependent electric field with period Tac- Since the coefficients CnK{i) in Eq.(p^) 
are periodic in t and A{t + Tac) = A{t) + A(Tac), we find that 
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(22) 



where AX = eA{Tac) /Tt-. The quasienergy e for a given solution Cxjit) to Eq.(Q) is only given modulo huJac, where 
ujac = '2.1^ /Tac- The substitution e ^ e + huaJ , I = integer, results in CKj{t) CKj{t)e^'^'"'* . Thus the Floquet state 
changes only by a time-dependent phase factor. 

Let us now consider two special cases for AK. First, if AK = modulo 2Tr/d, as for a harmonic field, the condition 
Eq.(p2|) is equivalent to Eq.(p7|) in the static case. All momenta in the BZ are then independent. 

A special case arises if AK = ^ | with p and q integers . Now a finite number of crystal momenta couple and 
the quasienergies are determined from 

CnK{t + qTac) = aK(i)e-"«^"=/'' . (23) 

The quasienergies are now only defined modulo hujac/q- The Brillouin zone collapses to [— Tr/qd, Tr/qd] and the 
quasienergy spectrum consists of fractional quasienergy ladders. 

We have so far assumed that the system stays in equilibrium until t = 0, when it is instantaneously coupled to the 
electric field. This is obviously an idealization, and it is natural to ask whether a finite switch-on period will change 
the general results. We have repeated the above analysis, but allowed for a finite switch-on period [0, To]. The upshot 
is that all formal results for the quasienergy spectrum still apply, however, one must change the initial momentum 
label K. The new momentum values, K ^ K — eA{To)/h, correspond to a translation in reciprocal space [ p7| . 

B. Tight-binding model 

In the rest of this paper we consider a superlattice model with only two bands, n — a,b. The Hamiltonian is 
= ^{ [A[5 - F{t)X'' - F{t)ld] ajai + [Ag - F{t)X'' - F{t)ld] h\hi 

-^{a\^^ai + 4ai+i) + ^{b\^^bi + b\bi+i) - F{t)X-\a\bi + fota,)} • (24) 

The terms involving the constants and X^ describe a site-independent shift, due to the electric field, of the 
unperturbed energies Aq and Aq, while the band coupling is described by the term involving the constant X°-^. The 
remaining terms in the Hamiltonian involve the site-dependent potential energy —F(t)ld and the hopping between 
nearest-neighbor sites. Due to the terms involving AT" and X^ the model represents a generalization of the two-band 
model which was first solved in [ p8| . The reason for introducing these parameters is that we can construct, as shown 
below, a mapping between the tight-binding model and a general two-band model (when the band couplings are 
described with the parameters i?", i?^, and i?"*"). The various energy parameters describing the two-band superlattice 
are summarized in Fig. |l|. 

We shall first consider static electric fields. Then the Schrodinger equation, -ff^u ~ eu, in the site-representation 
is solved by making the ansatz u — '^liu'^aj + u^6|)|0), and projecting out the /:th component, which results in the 
following coupled equations for the expansion coefficients: 

(Ag _ FX'^ _ Fld)ut - ^{uU + <+i) - FX^'ul = 6< 

(Ag - FX" - Fld)u'i + ^{ul, + u1+,) - FX'^'ut = euf . (25) 
For vanishing field the spectrum consists of two bands 

e^'^K) = AS''' T (A°'V2) cos(Xd) , (26) 

corresponding to the plane- wave solutions uf ~ e*^''^ , u'^ — and u° = , = giKid^ respectively. When the field 
is finite, we depart from the procedure used in [E8l by applying the discrete Fourier transformation 



K 



A- = l^E<^"^™' (27) 
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directly on Eq.(p5|). With the notation u'^ = Pj^e"^^/^ the result is 

e^iK) - FX'' -FX'''' 
-FX"'' e''{K) - FX 





(28) 



and 



PK+2./d - i'/^e-"^^/' . (29) 



The form of Eq.(p9[) is analogous to Eg. (po[) . Furthermore, Eq.(28) in the tight-binding model has the form of Eq.([l9|) 
if the bands are identical. In addition, the band couplings in the general model must be time-independent and satisfy 
X" ^ R"-, X'' = B!' and X"*" = R"^ = R''". The last equation requires that X"'' is real. If these conditions are fulfilled, 
the tight-binding model describes exactly the energy spectrum Eq.(^l|). The condition X"" — i?" = X'' — i?^ = can 
be relaxed in the case where only the difference spectrum between the two WS ladders is required to be exact. Then 
it is sufficient that X'' - R'' = X'' - R'' = AR g]. 

We now return to the solution of Eq.(p8|). The diagonal matrix elements in this equation can be eliminated by 

„__.„.__o . ^ . ^„ i'^'^-P^ndK ^ Q^^^ ^^^^^ ^^^^ g pa 



- I 



dK^ F dK 



(X"'') P^- = , (30) 



where we defined ujij{K) = uj^{K) + FX_/fi. Here the band separation is given by huo^{K) — e''(if ) — e°(/ir), and we 
also introduced X- = X" — X'' . The boundary condition is 

Let us next consider a part / of the Brillouin zone, where uJd{K) is a weak function of K (i.e., the bands are almost 
flat). Specifically, we demand that Aujd{I) = Jnax[u!d(K)] — mm[LLjd{K)], K G I, is small compared to (w^)/ + FX^/h 
(the average is calculated over /). In this case Eq.(|30|) is easily solved, and one finds 

Here the frequency uji is given by ujf — uj'^ + uj'^ with ujc — 2(|X'^''|/(i)ijjB. Let us define the mean band separation 
loq = {u>d)BZ and cDo = "^o + FX^/fi. The integration interval can be extended to the whole Brillouin zone when 
AuJd{BZ) ^ cDo- The solutions Eq.(^2|) then correspond to the spectrum 

+ A?, + Ak X+ hQi 

e± = " ^ " - F^ ± -^+phu;B , P integer , 

nf=ujl+^l, (33) 

where X^ — X°- + X''. Thus the two interpenetrating WS ladders are positioned symmetrically around the band 
midpoint, translated by — FX+/2, and the ladder separation is hfli. The exact form of the spectrum (^) can also 
be obtained in a much more direct way by looking for localized solutions to Eq.(p5|). We restrict the solution to be 
confined within the p:th unit cell. Solving 



Ag - F{X'' +pd)-e -FX"'' 

-FX'''' Af5 - F{X'' +pd)-e 



(34) 



we find the spectrum (|33|). The difference spectrum obtained under the flat band (or localization) condition will 
reappear in the following section when we discuss the electronic response in a static field in the case of nearly flat 
bands. 

Next let us turn to time-periodic fields, when the right hand side of Eq.(p5|) acquires the additional terms ih-^u^^''\ 
Following Zhao et al. |^ we substitute 

u'{ = Q'l exp{i[ei + eA{t)ld]/h} , (35) 

and Fourier transform, whence 
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^From the Fourier transform of (^5|) we get the boundary condition 

QK{t + Tac) = Q''K+AKe-''^-'\ (37) 

with AK = eA{Tac)/h, as above. A comparison of Eq.(|6|) with Eq.(|l|), and Eq.(|^ with Eq.(|2|), shows that 
in order to obtain a correspondence we must set the Fourier variable K equal to the initial crystal momentum in 
the general model. The (difference) quasienergy spectrum in the tight-binding model and the (difference) spectrum 
obtained in the general model are then seen to be identical under the same conditions as in the static case. 




III. RESPONSE 



We have thus far examined the general properties of the energy spectrum and the wave functions of a superlattice 
in an external field. The next step is to compute the physical quantities, such as densities or drift velocities. To this 
end we must construct an appropriate kinetic equation. The simplest approach would be to consider the Boltzmann 
equation. However, since we are interested in the properties of systems with two or more occupied minibands, where 
a coherent band-to-band transfer plays an essential role, it is necessary to use a method, which allows one to consider 
quantum mechanical tunneling. In the following sections we shall use the density matrix description. 



A. The density-matrix equations of motion 

In order to describe the coherent transport we shall start from the equations of motions for the density matrix in 
the limit where collisions may be neglected. Such equations of motions were derived by Krieger and lafrate p5| , p6| . 
Thus we take as our starting point Eq.(51) of Ref. adapted to the case of two minibands. The density matrix is 
defined as 

PnKn'K'it) = {^priK [x , t)\p{x , t)\i>n' K' {x , t)) , (38) 

and for coherent motion its diagonal elements (p^ = PaKaK etc.) satisfy 

= \Kei^h{K,t) j\*{K,t') [p'K{t')-pW)] rfij (39) 

dt 2' 

where 



Re h*{K,t) / h{K,t') [p'i,[t') - p^(t')] dt' , (40) 



^ 2F{t)R->'[k{t)] f^\^,[kit')]+F(t')R^lkit')]/n]dt' 

' h 



(41) 



and R-{k) = i?''(fc) - R''{k). Defining p±{K,t) = pj^{t) ± p]^(t) we get from Eq.(|§,|^ 

p+{K,t)^Q (42) 
P-{K,t) ^ -Re\h{K,t) h*{K,t')p^{K,t')dt'^ , 



(43) 



where we defined p = dp/dt. Eq.(^2|) is simply a statement of particle number conservation, while the integro- 
differential equation (43) determines the kinetic properties of our system. 

As it stands, Eq.(4^) is unsuited for analytical applications, because of the non-Markovian "collision term" 
(presently, the "collisions" consist of band-to-band transfers, dressed with the effects of the external fields). However, 
it can be solved numerically by the technique described in Ref. We now prove that Eq.(^3|) can be reduced 

to a third order differential equation, and then handled by standard numerical techniques, if we assume that u 7^ 
except in equilibrium. In addition it will be shown that the form of this equation is invariant after a finite switch-on 
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period [0, Tq]. At the end of this time interval the response belonging to an initial momentum K depends only on the 
position k{To) of the semiclassical motion of and on the boundary conditions of p- at t = Tq. 

Since the different ivT-components do not mix in (^), we can keep the initial momentum K fixed. In order to ease 
the notation we write Eq.(|4l|) as ft. = ue**^. Then the relation uh = (ii + iu(f))h holds. Using this, and differentiating 
Eq.dl) we get 



t 

3 , 



up- = iip- — utpHc <^ih J h*p- dt' j — u p- . (44) 

An auxiliary function w is next defined by w = — Im It satisfies the equation 

uw — iiw + u(f)p-. (45) 
We proceed by expressing w in terms of another auxiliary function, w = uw. Using uw = 4>p-, we obtain 

(46) 

This allows one to rewrite Eq.(^) as 

p- - -p- + u^p- +u<p [ — dt' = 0. (47) 

u Jo y- 

Finally, dividing Eq.(^) by ucj), differentiating and multiplying by u^cjP', we arrive at the third-order differential 
equation 

P3P- + P2P- + PiP- + Pop- = , (48) 

where 

P3 = u^(t> 

P2 = —u{2uU(j) + U(j)) 

Pi = {—uil + 2v? + u'' + M^0^)0 + uu(^ 

Pa = u^{u<j) - u(f)) . (49) 

/,From Eq.(p3|) and Eq.(p7|) above we get the boundary conditions p-{K,0) = and p-{K,0) — —v?p-{K,Q). These 
equations are readily numerically integrated, and a few special cases were described in our recent Letter pi|. 



1. The current 



The results derived above can be used to obtain an expression for the drift velocity v{t). The result can be expressed 
as a sum of a semiclassical, intraband term, and a term with quantum mechanical origin, which incorporates Zener 
tunneling. 

Omitting the details, we state our result: 

v{t) =Y.PK{i)^Kit) + 2^ [RepJ^nO Re<(0 - ln^p^^{t) lm<(0] , (50) 

vK K 

where 

1 de'^[k{t)] 

<(t) - -iujd[m]R''\k{t)] . (52) 
The non-diagonal density matrix elements can expressed in terms of the diagonal elements, 
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u(t)p^^{t) = -^-h{K,t) I h*{K,t')p_{K,t')dt' . (53) 



2 

Using the auxiliary function w defined in the previous section, we obtain 



Re{uit)p''^{t)} / ^dt' (54) 

^ Jo " 

Im {u{t)p''^ {t}} ^-^p_{K,t). (55) 



If R is real, Eg. (|50|) can be further reduced, and one finds 

uK ^ ^ K 

In Eq.(|5^) or (^6|) the first sum is the semiclassical contribution to the drift velocity, while the second sum accounts 
for Zener tunneling. Since p^(i) + Pxit) = const, (determined by the initial conditions), we can express Eq.(^ in 
terms of p-(K,t) alone. Thus, a numerical solution of the third-order differential equation derived in the previous 
section, Eq.(Q), is sufficient to determine the full time-dependence of v{t). This, however, must be done for each 
initially occupied if-value, and we leave applications for future work. 



B. Zener tunneling 



The numerical calculations of Zener resonances based on ( [48D were discussed in our recent Letter Since 
tunneling between bands plays an important role in the interpretation of our results we now consider the condition 
for Zener tunnehng in the present two-band model. The electrical field is assumed to be static F{t) = F{> 0). 



1. Analytical considerations 

a. Weak coupling limit. We shall show how a two-band version of Eq.(|l^) can give Zener resonances, following the 
perturbative method used by Mullen et al., Ref. [Q, which is valid for weak band coupling. The two-band system is 

(57) 

The system is linear and time-periodic with period Tb- The solution {C'^{NTb),C'1^{NTb)Y can then be obtained 
from the initial state (C^(0), C^(0))* by applying the time propagator Uq N times on the initial state. The propagator 

Uq is found as follows: By defining C^{t) — C^{t)e^'^ -fo ''^ the diagonal in Eq.(^ is ehminated, 

(58) 







where G"'' — i y: e jo . We work in the limit of small band couplings, \G°'^{t)TB\ ^ 1. To first order in 

r 



= J^"" G°-''dt' the initial state (1,0)* develops in a Bloch period Tb to (1, -(r°^)*)*, and the orthogonal state 



''dt 

(0, 1) to (r°'', 1)*. The time propagator is then 

Uo = e-r(--^-'--«'»-'[ (59) 

\ -(r"''Af)* M* I 

..^^^^^ ..^ — ^ '^'^ . According to ref. the resonance condition is that there should be a phase difference 27rp, 
p integer, between the diagonal elements. In our case, the resonance condition is therefore jj"^ Cjd dt' — 2'Kp or 

(60) 
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b. Static flat-band limit. Let us now discuss the case where the band couphng is not assumed to be smalL For a 
static field Eq.(|48l) reduces to |^ 



(61) 



We restrict ourselves to the case when tui — lu^i — ujq ujq. It is furthermore assumed that the band couplings R are 
independent of time and we denote X" = , X^' ^ R'' , and X"'' = i?"*" = R''" . We have 



u = 2{\X"'''\ld)uB = uJc 

(j) = -(tJo + Wi) . 

Eq.(|l]) is then 

— (cDq + ^l)P- + ^1/9- — (i^O + ^l) i^c + ('^O + t^l)^) P- + ^l^lP- = 

We write Eq.(|6^) in the form 

(Lo + Li + Lr){pa + pi + Pr) = , 
where and , /i = 0, 1, r , are of the zeroth, first, and higher order in oji. We get 

Lap- = -^qP- - ^a^ip- 

LiP- = —LUip- + ioip- — uji{^f + 2ujI)p^ + uiuoIp- 



(62) 
(63) 
(64) 



(65) 



The boundary conditions are poiK,0) = p_(X,0), pi(X,0) = Pr{K,0) = 0, p^(X,0) = 0, and p^,{K,0) 
-u)1p^{K, 0). The zeroth order equation LoPo = is readily integrated: 



Po + ^fpo = t^oPo(O) . 



(66) 



Using the zeroth order equation and Eq.(66) in the first order equation L^pi + iipo = 0, and integrating, we find 



pi + Q-i pi 



-Wo 



uji (po - po(0)) + / wipo rfi' 



(67) 



It is now clear that po only contains the frequency fi; . If wi (/f) is an analytic function then the spectrum of oJi is 
puJB , P integer. From Eq.(p7[) we see that the Fourier spectrum of pi consists of the frequencies fi; -\-puiB , P integer. 

Let us turn to the particular case with cosine bands, Eq. (p6|) . As shown in the previous section the cosine bands 
can be generated by a simple tight-binding model. For finite field values the energy spectrum is given by Eq.(^3|). 
Performing the integrals we obtain 



Po{K,t) = 
Pi{K,t) = 



Po{K,0) 



;^ cos flit) 



woiA-,+A\)uj^^Po{K,0) 



^ujb cos Kd + ujb cos{ujBt + Kd) + {fli — ujb) cos Kdcos flit — fli cos{flit + Kd) 

(68) 



COs[{fll+UJB)t + Kd] - ^' ^ C0S[(»; - UJB)t ~ Kd] | 



where fli was defined in (p3|). Thus to first order in uji the response contains the first intra WS transition lob and the 
inter WS transitions fli , fli^LOB- 

The structure of the perturbative solution ( p8| ) seems to suggest that the response frequency spectrum always 
corresponds to intra and inter transitions between (quasi) energy ladders. This expectation is proven analytically in 
Appendix A for static fields. The numerical examples in also support this conclusion for time-periodic fields. 
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2. A numerical example 



In our recent Letter pl[ | we presented examples of a numerical solution of the equation of motion for the density- 
matrix p^{K = 0,t). It was found that at certain, relatively sharp defined values of the applied field the time- 
dependence of /9_(0,t) exhibited a full inversion, i.e., its values ranged from -1-1 (the initial value corresponding to a 
carrier in band a) to -1 (carrier in band b). The sharp resonances, corresponding to complete band-to-band transfer, 
were termed "Zener resonances" . These features were shown to be intimately related to the Wannier-Stark ladder 
structure of the system: the resonances occurred when the level differences of the double ladder (corresponding 
to the two-band system) was at minimum. In the present section we describe the results of extensive simulations 
which generalize our earlier results in the following two respects: (i) We allow the interband coupling, described by a 
parameter f = X°'''/d, to vary from zero to large values; and (ii) We examine the equation of motion for 41 different 
iiT-points, covering the entire Brillouin zone. 

The results of our numerical work, which consists of three separate calculations, are summarized in Fig. ^. The 
first calculation consists of an evaluation of the local minima of the normalized WS ladder separation as a function 
of field strength. These are represented by diamonds in Fig. ^j. The normalized WS ladder separation is constructed 
in three steps: (i) Compute the differences of the energy levels in the two WS ladders; (ii) Divide the differences by 
Fd; and (iii) Select the part of the normahzed difference spectrum that belongs to [0,0.5]. It is seen from Fig. ||that 
by increasing the interband coupling the local minima are shifted monotonously towards higher field values and that 
they vanish in succession. The local minima are numbered from the high field end at zero band coupling. The first 
five local separation minima are depicted (when they exist) in Fig. H. 

The second part of Fig. |^ is a comparison of part one with the local minima of the WS ladder separation in the 
case where the double ladder is given by Eq. (|3^) . The normalized difference spectrum Cdif = ±(e~'' — e~)/Fd is 



€d^f = g ± W 4^2 + ( — j ; q integer . (69) 

In our numerical example we have put X"" — = 0, thus X^ = 0. The normalized WS ladder separation have local 
minima (equal to zero) whenever q — ^4^^ + {A'^'' / Fd)^ = or 

1 1 



Fd M^^^^- 

In the limit of zero band coupling we get 1/Fd — g/A"'' in exact agreement with the correct result. In the high 
field limit the local minima should be at ^ = q/2. As seen from Fig. ||, the flat band (or localized) model tends to 
overshoot the exact findings, but it provides an important hint about where to find the local minima. In general the 
flat band model is not formally correct: the numerical calculations show that the two WS ladders never coincide (level 
repulsion) for non-zero band coupling. But even in our case, where the bands are far from being flat, it estimates well 
the local minima. 

The last part of Fig. ^, denoted by error bars, is a calculation of Zener resonances, based on examining the time- 
development of {K, t) . As showed in our Letter |^l[ the number of oscillations in p_ {K, t) within a plateau increases 
when the field is decreased. A plateau is better defined at the centre than at the terminal points where the deviation 
from a constant value of p-{K, t) is largest. Thus for Zener resonances with high index it is not sufficient to calculate 
^^"■(if). Instead we define the centre of a plateau as the time where the oscillation amplitude of p-{K,t) is at 
local minimum. A Zener resonance is then a field value at which any of the plateau centre values have reached — 1. 
For a fixed interband coupling the first five Zener resonances have been estimated (when they exist) for 41 evenly 
distributed K points in the Brillouin zone. Some dispersion was detected, and this is indicated by the height of the 
error bars in Fig. |2 In the weak coupling limit we find zero dispersion and the Zener resonances coincide with Eq. (|60|) . 
As seen from Fig. A there is an excellent overall correspondence between local minima of the normalized WS ladder 



separation on one 



land, and the Zener resonances on the other hand. 



IV. CONCLUSION 



Our main results can be summarized as follows, (i) The wave function for an electron in a periodic potential, such 
as in a semiconductor superlattice, under the influence of a (time-dependent) uniform electric field is expanded in 
terms of accelerated Bloch states. The properties of the expansion coefficients allow us to give a unified treatment 
of several different results, obtained by a variety of techniques in the literature, for the superlattice. In static fields. 
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the existence of the Wannier-Stark ladder can be deduced, and in the time-dependent case, with a periodic external 
field, a fractional quasienergy spectrum emerges for certain values of the external field, (ii) A mapping is constructed 
between a two- band tight-binding model, and the general system considered in (i). This allows one to make several 
general statements concerning the properties of the two-band system, both in static and time-dependent situations, 
(iii) A density-matrix equation of motion is developed, and the resulting third-order differential equation is shown 
to be well-suited for both analytic and numerical studies of Zener tunneling in the two-miniband system. We also 
give new numerical results for Zener resonances, identified in our earlier work and interpret them analytically 
applying the tools developed in this work. 



APPENDIX A: THE RESPONSE FREQUENCY SPECTRUM 



We shall demonstrate that the response frequency spectrum in a dc field consists of elements corresponding to intra- 
and inter- WS transitions. In Sect.|| we found the linearly independent solutions 

^Kj{x,t) = e-"^^*/'"'^cj^,„0„K(a:) (Al) 

n 

to the Schrodinger equation by using the A^-dimensional Floquct theorem. Now expanding a general solution as 

i>{x,t) = ^ aKji>Kj ix,t) , (A2) 

we obtain for the diagonal element PnK{t) of the density operator p{x,t) — \tli{x,t)){il}(x^t)\: 

PnK{t) = {lpnK{x,t)\p{x,t)\lj}nK{x,t)) 

= aif'j'aK"/'CK'/„'c;^»jv,„,,e-'('^'^'-'^''^'')*/'"'(V'„K ■ (A3) 

K"j"n" 

The ABS are time-periodic functions with period Tg. Thus PnK{t) contains only frequencies which are a subset of 
energy differences between the WS ladders. 

In Section III the gauge was changed. We now show that the response frequency spectrum still consists of energy 
differences between the WS ladders for the scalar-potential Hamiltonian. In the new gauge the Schrodinger equation 
is on the form 

HA{x,t)^{x,t)^zh^^i^^ . (A4) 
We expand the wave function directly in the ABS %l){x,t) = J2nK -^riK{t)ipnK{x,t) and find the system 

^^ dB,,K{t) ^ _ F{t)Y,Rnn'{k{t))B^>K{t) , (A5) 

n' 

which is identical to Eq.(p^. The integration method in Section II gives us linearly independent solutions 

^kM^ '^Kjn{t)i^nK{x, t) . (A6) 

n 

Expanding the wavefunction ip{x^ t) — X^Kj ^Kj^^xji^^ 

PnKit) = <^Kra*^^.,e^^--^—^"'>'P%^,^{t)P%^„:{t) . (A7) 
j'j" 



The diagonal elements have a different appearance as compared to Eq.( A3), but the conclusion concerning the response 
frequency spectrum is the same as before the gauge change. 
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Ae(K) 




FIG. 1. Definitions of various energy parameters for tlie two-band superlattice. Aq'*" are tiie band mid-points, 
A"'' their widtlis, and Tiu)d{K) = ^b{K) — (-a{K) is the band separation at K. Also in the text are used: (i) 
AwdC/) = m&-K[ujd{K)] - mm[LJd{K)], K £ I; (ii) u>d{K) = uJd{K) + FX-/h (here X- = X" - X'', where X' are defined 
in Eq.(p^); (iii) {uid)i, which is the average of uJd{K), over a given part of the Brillouin zone /; (iv) = {uJd{K)) bz , i.e., the 
average band separation over the entire Brillouin zone, and, finally (v) uio ~ + FX-/h. 
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FIG. 2. The first five Zener resonances as a function of the band coupling parameter ^ = X"'^ /d. The diamonds represent 

the field values at which the distance between WS ladders has a local minimum. The solid lines give the local minima for the 
flat band model. The 'error bars' indicate at which field strengths the electronic response shows a Zener resonance. The finite 
widths are a consequence of dispersion. See the text for further discussion. The superlattice parameters are X- = 0, d = lOnm, 
AJ = 14meV, A\ = 14meV, and A"^ = A^ - Ag = 20meV. 
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